DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Sediment  Transport  at  Density  Fronts  in  Shallow  Water: 
a  Continuation  of  N00014-08-1-0846 

David  K.  Ralston 

Applied  Ocean  Physics  and  Engineering,  MS  #12 
Woods  Hole,  MA  02543 

phone:  (508)  289-2587  fax:(508)457-2194  email:  dralston@whoi.edu 

Award  Number:  N00014-13-1-0197 


LONG-TERM  GOALS 

The  goal  of  the  overall  research  is  to  quantify  through  observations  and  modeling  how  density  fronts  in 
shallow  estuarine  flows  impact  the  mobilization,  redistribution,  trapping,  and  deposition  of  suspended 
sediment.  This  current  award  is  a  continuation  of  the  original  YIP  award  in  order  to  deliver  a 
remaining  increment  of  funds  that  had  not  been  allocated.  The  focus  of  this  continuation  is  to  examine 
how  stratification  on  the  tidal  flats  varies  in  parameter  space  as  a  function  of  dimensionless  numbers. 

OBJECTIVES 

The  objectives  of  this  continuation  award  are  to 

analyze  results  from  a  high-resolution,  3 -dimensional,  finite-volume  hydrodynamic  model  of 
the  Skagit  tidal  flats  to  examine  processes  leading  to  the  creation  and  destruction  of 
stratification, 

-  use  the  Simpson  number  to  represent  the  balance  between  tidal  straining  and  tidal  mixing  to 
place  the  physical  processes  on  the  tidal  flats  in  a  more  general  context  of  estuarine  and  coastal 
stratification  dynamics 

APPROACH 

The  research  builds  off  an  extensive  observation  and  modeling  field  program  on  the  tidal  flats  of  the 
Skagit  River  delta.  In  the  field,  we  measured  velocity,  stratification,  and  suspended  sediment  at  high 
resolution  in  shallow  flows,  tracking  the  evolution  of  salinity  fronts  through  the  tidal  cycle.  A  major 
field  effort  occurred  in  June  2009  on  the  Skagit  Tidal  flats  in  Puget  Sound,  coordinated  with  other 
researchers  in  the  Tidal  Flats  DRI.  Focused  observations  of  the  shallow  density  front  and  its  evolution 
through  the  tidal  cycle  were  complemented  by  a  large  scale  array  of  moored  instruments  deployed 
during  the  same  period  (along  with  Geyer  and  Traykovski).  These  observations,  along  with 
observations  collected  by  other  investigators  (Raubenheimer  and  Elgar)  on  other  parts  of  the  Skagit 
tidal  flats  are  used  to  assess  the  validity  of  the  numerical  model  results. 

In  parallel  with  the  observations,  we  developed  a  numerical  model  of  the  Skagit  tidal  flats.  The  model 
uses  the  Finite  Volume  Coastal  Ocean  Model  (FVCOM),  but  was  modified  to  incorporate  recent 
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advancements  in  sediment  transport  modeling  with  the  Community  Sediment  Transport  Model  System 
(CSTMS).  The  unstructured  grid  of  FVCOM  allows  the  model  to  simulate  conditions  the  Skagit  flats 
with  enhanced  grid  resolution  near  the  observations.  The  observations  were  used  to  calibrate  the  model 
and  to  evaluate  how  well  the  model  represents  sharp  salinity  gradients  at  fronts,  both  across  the  tidal 
flats  and  at  lateral  fronts  coinciding  with  channel-shoal  bathymetry.  Collectively,  analyses  of  the 
observations  and  model  were  used  to  quantify  how  local  frontal  processes  on  scales  of  10’s  to  100’s  of 
meters  impact  retention,  redistribution,  and  export  of  sediment  over  tidal  flats  on  scales  of  kilometers. 
Model  results  from  simulations  of  the  observation  periods  in  2009  and  idealized  tidal  and  river  forcing 
conditions  provide  the  basis  for  the  analysis  of  stratification  on  the  tidal  flats  in  this  current  project. 

The  analysis  will  be  a  continuation  of  a  collaboration  with  V.  Pavel,  a  recent  graduate  of  the  WHOI 
joint  program,  and  her  advisor  B.  Raubenheimer. 
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Figure  1.  Bathymetry  of  the  Skagit  Bay  tidal  flats  (left)  with  a  zoom  on  the  study  area  on  the 
southern  flats  (right).  Red  dots  indicate  frame  locations  and  lines  show  across-flats  and  across- 

channel  survey  lines. 


WORK  COMPLETED 

Extensive  field  efforts  and  model  development  occurred  under  related  projects,  but  here  the  focus  has 
been  on  assessing  the  proceses  governing  stratification  on  the  tidal  flats  (Fig.  2).  Data  collected  during 
field  campaigns  on  the  Skagit  tidal  flats  in  2009  were  compared  with  the  model  results  for  those 
periods.  Several  manuscripts  have  been  published  recently  related  to  this  effort.  Topics  included 
quantifying  impacts  of  estuarine  and  fluvial  process  on  sediment  fluxes  (Ralston  et  al.,  2013),  an 
analyses  of  generation  mechanisms  for  tidal  asymmetries  on  the  tidal  flats  (Nidzieko  and  Ralston, 
2012),  and  an  evaluation  of  effects  of  complex  topography  on  wind  correlation  length  scales  and 
implications  for  coastal  ocean  modeling  (Raubenheimer  et  al.,  2013). 

Recently,  we  have  worked  on  applying  to  the  model  results  some  analytical  approaches  to  assess 
mechancisms  of  stratification  creation  and  destruction  on  the  tidal  flats.  This  work  is  in  collaboration 
with  V.  Pavel,  who  recently  completed  her  Ph.D.  at  WHOI  working  with  B.  Raubenheimer  on 
stratification  dynamics  on  the  Skagit  tidal  flats  based  predominantly  on  field  obsercvations.  The  model 
provides  the  opportunity  to  generalize  her  observational  results  to  other  parts  of  the  tidal  flats  and  to 
link  the  Skagit  findings  to  coastal  settings  with  similar  forcing  regimes. 
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Figure  2.  Across-shore  transects  of  salinity  on  the  Skagit  flats,  from  observations  (left)  and  model 
(right).  Top  2  rows  are  during  a  flood  tide,  with  a  relatively  strong  horizontal  salinity  gradient  and 
little  stratification  in  the  front;  in  both  the  model  and  observations  a  near  surface  layer  of fresher 
water  can  be  seen  behind  the  front.  Bottom  2  rows  are  during  an  ebb,  when  stratification  was  strong 
even  in  water  as  shallow  as  ~1  m.  Location  of  transect  is  shown  in  Fig.  1. 


RESULTS 


From  the  Skagit  model  results,  we  have  found  that  the  rate  of  change  of  stratification,  quantified  as  the 
integrated  potential  energy  anomaly  ®  (Simpson  et  al.  1990),  varied  with  the  Simpson  number, 


Si  = 


8j  dpi  2 

f  Po  J  dx H 


cdu3 


-,  where  g  is  gravity,  dp/dx  is  the  horizontal  density  gradient,  H  is  the  water  depth,  Cd 


is  a  drag  coefficient,  and  U  is  the  tidal  velocity  (Fig.  3).  The  Simpson  number  represents  a 
dimensionless  balance  between  the  creation  of  stratification  by  straining  of  the  horizontal  density 
gradient  versus  destruction  of  stratification  due  turbulent  mixing  associated  with  tidal  bottom  stresses 
(Stacey  et  al.  2001).  These  preliminary  results  were  from  a  narrow  set  of  conditions  focused  during  the 
period  of  observations  and  in  the  region  where  instruments  were  deployed.  The  on-going  analysis  is  to 
extend  these  results  over  a  longer  time  series  that  includes  a  wider  array  of  forcing  conditions,  and  to 
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assess  how  the  results  vary  spatially  with  mean  water  depth  and  proximity  to  the  river  mouth.  The 
Skagit  results  will  be  compared  with  results  from  models  of  other  estuaries,  including  the  Hudson, 
Merimack,  and  Connecticut  Rivers,  to  assess  how  the  stratification  and  the  Simpson  number  transitions 
found  on  the  tidal  flats  fits  into  estuarine  parameters  space.  If  the  Simpson  number  scaling  is 
appropriate,  then  the  dimensionless  scaling  should  hold  despite  the  differences  in  flow  environments. 
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Figure  3.  Rate  of  change  of  stratification  (d&/dt)  on  the  Skagit  tidal  flats  vs.  Simpson  number,  from 
numerical  model  results  (Pavel,  2012).  Type  1  tides  are  weaker  and  predominantly  semi-diurnal, 
while  Type  2  tides  are  more  energetic  and  have  a  greater  diurnal  component. 


IMPACT/APPLICATIONS 

Results  from  this  project  may  be  used  to  enhance  hydrodynamic  and  morphological  models  of 
estuaries  and  deltas,  with  applications  to  environmental  assessment  for  the  Navy.  Trapping  and 
deposition  of  sediment  associated  with  stratification  and  density  fronts  could  introduce  spatial  and 
temporal  variability  in  bed  consolidation  and  bathymetric  relief  on  tidal  flats.  The  project  will  also  help 
to  evaluate  the  skill  of  coastal  hydrodynamic  models  at  resolving  density  fronts  and  stratification, 
including  the  surface  expressions  that  can  be  assessed  with  remote  sensing. 

RELATED  PROJECTS 

The  work  here  is  a  continuation  of  the  YIP  awarded  to  Ralston  (N000 14-08- 1-0846),  and  is  closely 
linked  to  the  Tidal  Flats  DRI.  The  field  efforts  on  the  Skagit  were  done  in  conjunction  with  Geyer  and 
Traykovski.  The  model  and  grid  development  has  been  in  collaboration  with  Geoff  Cowles. 
Collaborations  with  others  involved  in  the  DRI  include  Raubenheimer  and  Elgar  (for  bathymetry, 
observations  for  model  calibration,  and  use  of  the  model  to  interpret  water  column  and  wind 
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observations),  Lerczak  (model  simulations  at  seasonal  time  scales),  Signed  and  Sherwood  (CSTM 
implementation),  and  Thomson  and  Chickadel  (bathymetry). 
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